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Towards hydrodynamics without an entropy current 
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We present a generating functional which describes the equilibrium thermodynamic response of 
a relativistic system to external sources. A variational principle gives rise to constraints on the 
response parameters of relativistic hydrodynamics without making use of an entropy current. Our 
method reproduces and extends results available in the literature. It also provides a technique for 
efficiently computing n-point zero-frequency hydrodynamic correlation functions without the need 
to solve the equations of hydrodynamics. 



Introduction. — Hydrodynamics is a generic effective 
theory, valid on distance scales much longer than the 
typical mean free path, and applicable to many diverse 
physical theories at finite temperature [T] . The equations 
of hydrodynamics are characterized by several response 
parameters which need to be specified for each particu- 
lar system. These response parameters are usually con- 
strained by a set of equalities and inequalities vifhich are 
conventionally determined by requiring the existence of 
a local entropy current with positive semi-definite diver- 
gence. In this Letter, we will systematically derive the 
equality-type constraints on the response parameters us- 
ing a variational principle. 

In the hydrodynamic regime, a relativistic system can 
be described in terms of a velocity field u^, normalized 
such that Ufj,u^ = — 1, and a temperature T. When there 
is a conserved U{1) charge, the corresponding chemical 
potential ^ provides an additional hydrodynamic degree 
of freedom. If the U{1) symmetry is spontaneously bro- 
ken, the emerging Goldstone boson also turns into a 
hydrodynamic degree of freedom. 

The energy-momentum tensor T^'' and the (non- 
anomalous) charge current J'^ may be expressed through 
constitutive relations in terms of the hydrodynamic vari- 
ables and their gradients. The kinematic equations for 
hydrodynamics then amount to energy-momentum and 
charge conservation, 



Df.Ti"' 



F'^PJ,, 



D^J>' 



0. 



.p, ^p. -v.. (1) 

In (IT]) F^-'^ is the field strength of a background gauge 
field A^ conjugate to J^. The covariant derivative D^ 
depends on a background metric g^^. 

In d spacetime dimensions, we decompose the energy- 
momentum tensor and U{1) current into scalars, vectors 
and tensors of the S0{d—1) C 50(^—1,1) symmetry 
preserved by u^. 






(2) 



where A'^" = g''" 



and u^r^" 



''u'^ is a projection matrix, q^'u^ — 
g^^T'^'' = 0. The scalars £", V, 



and Af along with the vectors q^ , i^'' and the tensor r^"^ 
may be written as local functions of the hydrodynamic 
variables and their derivatives. In the hydrodynamic ap- 
proximation, the constitutive relations (pi) can be written 
in a derivative expansion 0. 

Several considerations come into play in determining 
which tensor structures can contribute to the quantities 
in (pi). First, we note that there is an inherent ambiguity 
in defining the velocity field, temperature, and chemical 
potential out of equilibrium. We may always redefine 
T ^ T + ST, fi ^ fi + 6fi, and m^ ^ u'' -H (Su'' such 
that ST, 6fj,, and Su^ vanish in the absence of gradients. 
Such field redefinitions are called changes of frame. A 
canonical choice of frame is the Landau frame in which 
qP^ — 0, £ — e, and Af = p with e and p the energy and 
charge densities in the absence of gradients. Even after 
choosing a frame not all tensor structures are allowed; as 
it turns out, the existence of an entropy current together 
with the Onsager relations leads to restrictions on the 
allowed tensor structures [I]. 

The restrictions imposed by the existence of an entropy 
current are either inequalities or equalities. For example, 
in the presence of an electric field, the tensor decompo- 
sition of the current, J* = kA'^'S^^ + aE^, is subject 
to an equality-type relation, k — —aT. In this Letter, 
we systematically show how relations of this type are en- 
forced by the equilibrium properties of the theory in the 
presence of external sources A^ and g^^, and follow from 
a variational principle. As a result, we learn that the hy- 
drodynamic constitutive relations are constrained both 
by symmetry and by the need to consistently describe 
static equilibria with external sources. 

Equilibria. — A time-independent equilibrium configu- 
ration can be characterized by a constant timelike vector 
V^^, where V^ = (1,0) in suitable coordinates. Starting 
from a source-free equilibrium configuration we assume 
that finite sources can be turned on adiabatically while 
maintaining equilibrium. In other words, we will only 
be considering configurations in which the Lie derivative 
with respect to V^, £y, vanishes when acting on ther- 
modynamic quantities or sources. Furthermore, we will 



be studying configurations with a finite static correla- 
tion length. Thus, Euclidean correlation functions fall 
off exponentially at large distances, implying that zero- 
frequency Fourier-space correlators are analytic at low 
momentum. 

Correlators in the equilibrium configuration can be ob- 
tained by differentiating a generating functional with re- 
spect to the sources. Indeed, consider the set of all zero- 
frequency correlation functions, expanded to m^^ order 
in momenta about zero. We call these n-point func- 
tions truncated correlators. After a Fourier transform, we 
obtain approximate position-space correlation functions 
valid on length scales much larger than the correlation 
length of the system, much like a multipole approxima- 
tion characterizes a localized distribution on large scales. 
Integrating the truncated functions over sources leads to 
one-point functions, which will be local functions of the 
sources. These may be further integrated to obtain the 
equilibrium generating functional for truncated correla- 
tion functions 



W„,= 



d a::/![sources(x)]. 



(3) 



where C includes terms with up to m derivatives [3]. 

In order for W„i to be diffeomorphism and gauge in- 
variant, C must be constructed from local diffeomor- 
phism and gauge invariant scalars, possibly in combina- 
tion with V^. In addition, C can depend on observables 
that are local in space but non-local in Euclidean time 
such as the invariant length of the time circle in the Eu- 
clideanized theory L, and the Polyakov loops Pa of any 
U{1) gauge fields. Since £y = 0, we find L = ^\f—V^ 
and InP^ = fW^Ay,, where /? is the coordinate periodic- 
ity of the time circle [3] . We identify the temperature T, 
the chemical potential /i, and the velocity field u'' as 



r=l/L, ^t = lnF^/i, 






(4) 



The parameters T, /i, and u^ depend on position through 
A^ and F^. 

Suppose that there are iV„ scalar quantities at n*'' or- 
der in a derivative expansion. We will denote them by 
Sn.i, , s„^2j • ■ • J Sn,N„- For instaucc, in a theory contain- 
ing a single conserved current (corresponding to an un- 
broken symmetry) we have so.i = T and so,2 — M- The 
most general generating functional for truncated zero- 
frequency correlators is of the form 



W„ 



d'^x^ 



m N„ 



{so)Sr. 



= 1 i=l 



, (5) 



where the a„_i and P are functions of the zeroth or- 
der scalars which we denoted collectively by sq. In the 
source-free equilibrium state, all of the derivative contri- 
butions to (l5| vanish, so that Wm is the logarithm of the 



exact equilibrium partition function. Thus, we identify 
P with the pressure of the source-free equilibrium state. 
We obtain one-point functions of the energy- 
momentum tensor and conserved current by varying Wm 
with respect to the metric and gauge field. 



5A, 



(6) 



If we denote the set of n*^ order transverse vectors and 
transverse traceless tensors by d„ ^ and i„_i then, on com- 
paring (|6| to (pi), we find 



^n.i^n,i i 



n=0 i=l 

m Nn 

n=0 i=l 

m N„ 

-^ = EE^^ 



m Nr^ 






1=0 i=l 



n=0 1=1 
m N„ 

'?'' = EEt"'*^",^' 

n=0 i=l 

m N„ 

-Ml' _ \^ Y^ 5» it^" 



(7) 



1=0 1=1 



where the e's, tt's, 0's, 7's, 5's and 6''s are determined in 
terms of the a's and their derivatives. While the most 
general expression for the energy-momentum tensor and 
current takes the form ^, the fact that T'^'^ and J'^ 
have been obtained from a local generating functional 
implies that not all tensor, vector and scalar structures 
are allowed, and that relations of the form (l7| hold. 

The gauge and diffeomorphism invariance of Wm en- 
sures that the solution ^ provided by ^ will satisfy 
the hydrodynamic equations (fTj). Matching the thermo- 
dynamic theory to the effective hydrodynamic descrip- 
tion ([7]) gives the constitutive relations in a particular 
frame, which we call the thermodynamic frame. In this 
frame, the values for the temperature, chemical potential 
and velocity field remain unchanged from their equilib- 
rium definitions Q after the hydrodynamic equations 
have been solved. 

In what follows we will give several explicit examples 
of systems where the relations (It]) are obtained from the 
generating functional ^. Some of these systems have 
been analyzed in the literature by requiring the existence 
of an entropy current. In all our examples, the non- 
dissipative constraints obtained using the entropy current 
method match those obtained here. 

Example 1: Ideal superfluids. — We begin by construct- 
ing the generating functional and computing the one- 
point functions for a superfluid to zeroth order in deriva- 
tives (i.e. an ideal superfluid, see [SUB] for a brief review). 
In addition to the zeroth order scalars so,i — T and 
So, 2 = M '^e can, a priori, construct two scalars from the 
extra hydrodynamic degree of freedom associated with 
the Goldstone boson, ^^^^ = — ^^ and u^^^, where ^^ is 
the gauge invariant combination ^^ = — 9^(/' + A^. Since 
£y(/) = implies that u^S,p, = ^J,, only so,3 = ^^ is an in- 



dependent scalar. According to (Is]), the generating func- 
tional takes the form 



Wo 



-gP{T,ti,e 



Using ^ we find 



( J'^) = pu^^ - fe 



U^Cu 



(8) 



(9) 



with dP = sdT + pd^i + Ifd^^ imAe^Ts + pp-P T\. 
The expressions in (l9| precisely match those of an ideal 
superfluid in the notation of [5]. For ^ = we recover 
the standard expression for an ideal normal fluid. 

Example 2: Parity-violating fluids. — For parity- 
violating theories in 2-1-1 dimensions with a conserved 
U{\) charge, the zeroth order scalars are sq.i = T 
and So. 2 = M- At first order in the derivative expan- 
sion there are a priori three scalars, Df^u^, u'^d^T, and 
u'^d^fi. However, all three scalars vanish identically since 
Cv — 0. There are two non-vanishing pseudo-scalars, si 
and S2 defined in Table|l]which specify the magnetic field 
and vorticity. Thus, we have the generating functional 



Wi = 



d^3 



^ [P{T, p.) + aisi + a2S2] 



(10) 



The notation in ( 10 1 deviates slightly from (Is]) in that the 
coefficients of parity odd tensors are adorned with a tilde. 
Furthermore, since all the results in this example involve 
tensors with one derivative, we have also dropped the 
derivative index n. We use the same simplified notation 
for the quantities in ([t]). 

An exhaustive list of all possible first derivative tensors 
can be found in [5j. The constraint Cy — leads to the 
restricted list given in Table ll] Indeed, in equilibrium we 
find that 



d^T = -Ta^, df^p^ -paf^ + Ef_t, 



D„u„ 



-u^a^ +ujfj.^, 



(11) 



are satisfied identically, where we have defined a^ = 
u'^D^u'', and lu'"' = ^^^^^^{DpU„ ~ D„Up). Thus, the 
shear tensor a*^^ constructed from the projected, trace- 
less, symmetrized version of D^u^, vanishes as does the 



vector En 



-Tdp tj: . We note that both af"" and E^ - 



-Tdp^ 



contribute to dissipation, consistent with the claim that 
we are studying equilibrium states. 

By varying Wi with respect to the metric and gauge 
field and decomposing according to (pi) and (It]) we find 

TTi=j^=5i=0, 02=7i=a2-<5i , 0i=(5i=ai , (12) 

ei=T62=Ta'i+pai~ai, e2~T^2—Toi2+pa.2—'2.a2, 

where a prime denotes a derivative with respect to T and 
a dot a derivative with respect to p. 

An analysis of parity- violating hydrodynamics in 2 -I- 1 
dimensions based on a local version of the second law of 





1 


2 


pseudoscalars (si) 
vectors {vi) 
pseudovectors {vi) 


1 r^^^Pn, F 


^e^"'Pu^dpT 



TABLE I. Independent first-order data for 2-1-1 dimensional 
fluids. We have defined i5,i = F^^^u" . 
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3 


4 


scalars [si) 
vectors {vi) 
tensors {ti) 


R 

A^'^R^pu" 


u'^R^^u-' 

uj'"'a, 

-UpRf^^'^'^'^u^ 







TABLE II. Independent second order data. The expressions 
for a and uj are given by the inline expression following (fTTJ. 
R^vpa is the Riemann tensor and R the Ricci scalar. Triangu- 
lar brackets denote a projected traceless symmetrized tensor. 



{pu) 



2 ^fip^l^CT I 



A"" + A" 



^ff^'A^M' 



ct/S 



thermodynamics can be found in 9 . Those results were 
presented in the Landau frame with 



P = P -XBSi -XnS2 



XEV'^ 



XTV2 



XEv'i 



XTV2 



(13) 



where the ellipsis denotes tensors which vanish in the 
equilibrium states under consideration. Matching the 



thermodynamic result ( 12 ) to the Landau frame coef- 
ficients we find 



XB = TTl 



Xn = 7^2 



dP. 

Ik 

dP 



ei 



£2 



dP~ 

op 

dP~ 

*2 



Xe ^ Si- i?7i 



de dp 

XE= h- Rll , XT = ^2- Rl2 



XT =§2- Rl2 



(14) 



P), where ^ and ^ are evalu- 



with R = p/{e 
ated at fixed e and p respectively. From ( 14 ) we find 
Xe — XT = along with two relations among the four 
x's. These relations are identical to those found in [5] 
with /n — 0, M.B — <5i and M.n — 52- We emphasize 
that coefficients associated with tensor structures which 
vanish are undetermined by this method. 

Example 3: Second-order hydrodynamics. — In our fi- 
nal example, we consider a parity-preserving theory in 
d spacetime dimensions without conserved U{1) currents 
to second order in derivatives. There is only one zeroth 
order scalar, so,i = T. A computation similar to the one 
for the 2 -|- 1 dimensional fluid implies that there are no 
first order scalars and four second order scalars [10]; see 
Table [11} Using similar notation to the previous example, 
we will drop the derivative index from the quantities in 
^ and ([7|. 

Varying the generating functional W2 defined in ([5]) 
with respect to the metric and expanding the resulting 



energy-momentum tensor according to ^ and (l7| we 
find for d > 3, 

63 = T{a':j+a2—2a[) + 3{a2~a3)+2a4 , 

64 = T2(2a'/-4') + 2r(2a;-4)-ra4-a4 , 



(15) 



d-3 2(rf-2) ^ , 2 

a— 1 a— 1 d— 1 

d-5 2(d-2) , 
7r3 = (a3-a2)^ + ^33-rai, 

7r4 = ia,+T{a'2-2a[))^ ^^T^a'/, 

0i = -2ai, 03=4(a2-a3)-2Ta'i, 

6*2 = -2Ta[, 04 = 2(r2a"+r(2a'i-a2)-a4), 

71 = 2(ai+Q;2-a3), 72 = --2T{a[+a2-a'^) , 



where the tensor structures t^, Vi and Si are given in Table 
[TT| We refer the reader to [TT] for a comprehensive dis- 
cussion of second order tensor structures. In three space- 
time dimensions ( 15 ) still describes the expansion (l7| of 
the energy-momentum tensor, but the list of tensors is 
overcomplete. In particular, the tensors t^, and ti + ^2 
vanish so that only the combinations 9i — 62 and 64 ap- 
pear in T 



flU 



An analysis of the restrictions on response coefficients 
of 3-1-1 dimensional systems to second order in the deriva- 
tive expansion was carried out in [H] (see also [12] ) . The 
results were presented in the Landau frame where 



V =T (C2S1 + C3S2 + 6S3 + ^4S4) + • • • 

■^^- =T {n^e^" + K2tr + Aa^r + ^^^T) + 



(16) 



The frame transformation from the thermodynamic 
frame to the Landau frame is given by 6u^^ = q^ /{e -\- P) 
and 5T = (e — f )/e'. After carrying out this frame trans- 
formation, we find 



Tki = 01, Tk2 = 6*2, rA3 = 6*3, TA4 = 

dP ^^ dP 

oe oe 

dP ^^ dP 

oe oe 



(17) 



where ^ 



^^- Despite the fact that there are only 



four tti's, one can verify that there are five relations be- 
tween the eight coefficients above. These correspond pre- 
cisely to the five conditions on the response coefficients 
found in |llj . The remaining seven transport coefficients 
are undetermined either by requiring the existence of an 
entropy current or by the variational method described 
in this Letter. The results for d ^ 4 are new; the varia- 
tional method provides a simple alternative to the more 
onerous technique which uses the entropy current. 

Discussion. — In this work we have studied the impli- 
cations of the existence of an equilibrium state on the 



hydrodynamic constitutive relations. We have shown, us- 
ing three examples, how relations among response coeffi- 
cients which are canonically derived using a local version 
of the second law of thermodynamics emerge from prop- 
erties of the gauge- and diffeomorphism-invariant gener- 
ating functional ([s]). 

Varying the generating functional with respect to the 
sources leads to zero frequency n-point Euclidean corre- 
lation functions. Since the dependence of the generating 
functional on the sources is known explicitly, there is no 
need to solve the equations of hydrodynamics in order to 
compute the correlators. This significantly reduces the 
complexity of the computation. 

More importantly, the thermodynamic relations fol- 
lowing from the generating functional combined with the 
inequalities imposed on dynamical transport coefficients 
from positivity of spectral functions appear to repro- 
duce the entire suite of constraints implied by an en- 
tropy current with positive semi-definite divergence. Put 
differently, requiring the existence of a local entropy cur- 
rent with positive semi-definite divergence implies several 
constraints among the coefficients in the hydrodynamic 
equations. Some of these constraints appear in the form 
of equalities and others in the form of inequalities. As 
we have suggested in this Letter, the former constraints 
can be obtained by appealing to equilibrium thermody- 
namics. This suggestion may appear less surprising once 
we realize that equality-type constraints relate to dissi- 
pationless contributions to the energy-momentum ten- 
sor and charge current. It would be interesting to study 
whether all inequality-type constraints associated with 
dissipation follow in general from the positivity proper- 
ties of even n-point functions. 
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